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<D . Abstract 

fi: 

In (2+l)-dimensional hydro dynamic system with broken parity, there additionally exist 
the Hall viscosity and curl viscosity. The dual holographic model has been constructed by 
coupling a pseudo scalar to the gravitational Chern-Simons term in (3+l)-dimensional 
bulk gravity. In this paper, we investigate the non-relativist ic fluid with Hall viscosity 
and curl viscosity living on a finite radial cutoff surface in the bulk. Employing the 
non-relativistic hydrodynamic expansion method, we obtain the incompressible Navier- 
Stokes equations with Hall viscosity and curl viscosity. Unlike the shear viscosity, the 
ratio of the Hall viscosity over entropy density is found to be cutoff scale dependent, and 
it tends to zero when the cutoff surface approaches to the horizon of the background 
spacetime. 

PACS numbers: 11.25.Tq, 47.10.ad, 04.20.-q, 11.15.Yc 
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1 Introduction 

Over the past years there have been a lot of studies on the gauge/gravity dualities [1, 2, 3]. In 
the long-wavelength limit, the fluid/gravity correspondence relates the boundary hydrodynamic 
equations to the perturbation equations in the bulk gravity [4]. Recently, the (2+l)-dimensional 
system with broken parity has attracted much attention, and it has been holographically realized 
through coupling a pseudo scalar to the gravitational Chern-Simons term in the (3+l)-dimensional 
bulk gravity [5, 6, 7, 8, 9]. The dual gravity model is a generalization of the Chern-Simons modified 
gravity [10, 11], which has been widely studied in the context of cosmology, gravitational waves, 
and gravitational tests in the solar and earth systems, see e.g. [12, 13]. In the dual model, RR is 
the parity-violating Pontryagin density, which is a measure of the divergence of axial vector current 
of massless chiral fermions V^j 5 ^ in the gravitational anomaly [14]. 

The Hall viscosity T]a is a non-dissipative viscosity coefficient analogous to Hall conductivity 
which exists in the (2+l)-dimension fluid with time reversal symmetry breaking, it does not con- 
tribute to the entropy production of the fluid. In the quantum Hall fluids, the ordinary dissipative 
viscosities 77 and ( are absent at zero temperature, while the non-dissipative Hall viscosity r]A could 
exist [15]. It can be characterized as the rational numbers and provides a fundamental measure of 
incompressibility of the Hall fluid, and it is related to a Berry phase and proportional to the density 
of intrinsic angular momentum [16, 17, 18, 19, 20, 21]. It is also proposed that the Hall viscosity 
can be measured as the coefficient in front of q 2 term in the Hall conductivity in small wave number 
q limit [22]. 

In the non-relativistic limit of the (2+l)-dimensional hydrodynamics, the incompressible Navier- 
Stokes equations will be corrected by the parity odd viscosities, and it would be interesting to deduce 
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the equations of motion holographically. In this paper, we use the non-relativistic fluid expansion 
approach to study the (2+l)-dimensional non-relativistic hydrodynamics with Hall viscosity and 
curl viscosity, dual to gravity involving a pseudo scalar coupled to a topological gravitational Chern- 
Simmons term. By using the non-relativistic hydrodynamic expansion method associated with a 
flat cutoff surface, we formally solve the gravity equations and the pseudo scalar equation up to 
second order in the non-relativistic hydrodynamic expansion parameter. The Brown- York stress 
tensor on the cutoff surface is identified with the stress energy tensor of the dual fluid, from which 
we can read out the general analytic expression for the transports coefficients, such as shear viscosity 
77, Hall viscosity i]a and curl viscosity (a- The incompressible Navier-Stokes equations with Hall 
viscosity and curl viscosity are obtained from the gravity side. In addition, it is shown that the Hall 
viscosity over entropy density is cutoff scale dependent, and tends to zero when the cutoff surface 
approaches to the horizon of the background spacetime. 

This paper is organized as follows. In Section 2 the (2+l)-dimensional hydrodynamics with 
Hall viscosity and bulk viscosity in non-relativistic limit are investigated. Section 3 is the main 
part of this paper, where the incompressible Navier-Stokes equations with corrections are deduced 
from (3+l)-dimensional dual gravity. The case when the cutoff surface goes to the AdS boundary is 
discussed in Section 4. And some relevant further discussions are given in Section 5. In Appendix A, 
the holographic non-relativistic expansion procedure associated with a finite cutoff surface in (d+1)- 
dimensional gravity is briefly introduced. In this paper, we use the small letter = 1, 2, d — 1 
to denote the index of pure spatial coordinate x l , the Greek symbols /i, v, — 0, l,...,d — 1, to 
denote index of the ordinary space-time coordinate x^ ~ (r, x l ), and the capital letter M,N,= 
r, 0, 1, d — 1 to denote the index of bulk space-time coordinates x M ~ (r, x^ 1 ). 

2 Hydrodynamics with Hall Viscosity and Curl Viscosity 

In a (2+l)-dimensional parity violating hydrodynamic system, the energy-momentum tensor of the 
fluid with the first order gradient expansion can be written as [5, 6, 7, 8, 9] 

= eu"u v + pP" v - 2r)a^ - (QP^ - 2va<?T - (a^P^, (1) 

where the thermodynamic quantities e and p are the energy density and pressure respectively, which 
depend on the local temperature of the fluid and relate with each other through the equation of 
state p = p(e). The shear viscosity rj and bulk viscosity ( are the canonical transport coefficients, 
while the Hall viscosity i]a and curl viscosity (a arise from the parity violating effect. In a (2+1)- 
dimensional flat space-time background, the velocity = (1, /3 4 )/a/1 — (3 2 and the projection 
tensor P^ v = rf v + u^u v should be functions of the space-time coordinates x M = (r, x % ). The 
first-order gradient expansion tensors can be expressed as 

a*" = l -p^p^ (d aU(3 + d?u a - ep a ?) , e = v ^d,u u , (2) 

< = ^(e^W + e ^O> tt = -e^u^u p . (3) 

where we have used the convention e TXy = 1, and chosen the Landau frame T^ v u u = —eu^. 

In the thermal equilibrium state, we assume the fluid at rest has constant energy density Cq and 
pressure po- Through introducing the heat function per unit volume wq = eo +Po, we can define the 
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normalized fluctuations of the thermodynamic parameters as £ = (e — eo) / wq and V = (p — po)/wo- 
In the non-relativistic hydrodynamic limit, the following scaling appear [31, 34], 

PI = ePiiex', e 2 r), £ e = e 2 S(ex\ e 2 r), V e = e 2 V[ex\ e 2 r), (4) 

where e 1 is a small parameter. The velocity /?, could be regarded as a small vector fluctuation, 
while the thermodynamic quantities S and V basically come from the temperature fluctuation, and 
they would relate to each other through the equation of state p = p(e) of the fluid. Up to the second 
order of the non-relativistic expansion parameter e, one has 

u, ~ (-1, 0) + (0, A)e+i(/3 2 , 0)e 2 + O(e 3 ), (5) 

The non-zero components of the normal traceless shear viscosity tensor and the traceless Hall 
viscosity tensor are 

If (d X /3 x - dy(3y) (d X (3y + Oy (3 X ) \ O 3 

r A 1 ( (+d x l3 y + d y /3 x ) (—d x p x + dyPy) \ 2 , rnf ^ 



« 2 V (~dxl3 X + y f3y) (~d X f3y ~ dy(3 X ) ^ ^ ^' ^ 

the divergence and the curl of the velocity become 

6 = (d x p x + dy(3 y )e 2 + 0(e 3 ), Q = (d x f3 y - d y (3 x )e 2 + G(e 3 ). (9) 

The dynamical equations of the relativistic fluid are d^T^ = J 7 ", where J 71 " is an external force 
density, and we can define the normalized force density as f u = F v /wq. In the non-relativistic limit 
(4), d^T^ u = J 11 ' will reduce to the incompressible Navier-Stokes equations with Hall viscosity and 
curl viscosity, 

d l A = 0, d T fr + (3jd j (3i + d{P-u d 2 /3i - VAe ij d 2 Pj - Ue j %djP k = h, (10) 

where v = r]/wo, va = t/a/wo and £a = Ca/^o could be named as the kinematic shear viscosity, 
kinematic Hall viscosity and kinematic curl viscosity, respectively. 

3 Parity Breaking Hydrodynamics from Gravity 

In this section, we will deduce the (2+l)-dimensional non-relativistic parity breaking hydrodynamics 
described in the previous section from the (3+l)-dimensional gravity. 

3.1 The Dual Gravity Model 

The dual gravity model is described by the Einstein gravity involving a pseudo scalar coupled to 
the topological gravitational Chern-Simmons term. The action of the model can be expressed as [5] 

Smak = ! d A x^{R-2k + C cs )+ I d A x^£ e , (11) 



2k 4 
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where g = det qm n is the determinant of the metric, R is the Ricci scalar, A is the cosmological 
constant 1 , k 4 = 8tiGn with Gn the Newton's gravitational constant. The Lagrangian density Ccs 
and Lq are 

C cs = j9RR, C e = -\g MN V M 9V N 9 - V(9), (12) 

where A is coupling constant and 9 is the pseudo scalar field. The Pontryagin density RR is defined 

as 



pn_ r>M PQ r>N r>M PQ _ 1 PQAB nM 

nn — n N n MP q, n n — -e n tvab^ 



(13) 



where e MNAB is the four- dimensional Levi-Civita tensor in the bulk with the convention e rTXy = 
1 / y/—g. Varying the action with respect to the metric and the pseudo scalar respectively lead to 
the equations of motion 

Wmn = 0, Wmn = Rmn — ^9mnR + A^mjv + ACmjv — ^iT MN , (14) 

W e = 0, W e = 2k 4 (v 2 - ^p-^J ~ \RR, (15) 

where the stress energy tensor of the pseudo scalar is 

T e MN = -2^^ + g MN £g = V M 9V N 9-(^g AB V A 9V B 9 + V(9))g MN , (16) 

and the four-dimensional Cotton tensor Cmn is a symmetric traceless tensor of the second rank, 
which is defined through 2 [10, 11] 

C MN = -v s e SAB ^ M V A R^ +v S tR T{MN)S , (17) 
v s = V s 9, v ST = \7 s X7 T 9 = X7 {s \7 T) 9. (18) 

The Cotton tensor's divergence has a non-zero topological source as 

V M C MN = V RR= V d M J M , (19) 
o 4 

where J M origins from the 4-divergence of the gravitational Chern-Simons topological current 

jm = e ^ ( r^9prg A + ^r^rf 5 r^). (20) 

The gravity equations Wmn = simultaneously imply V m Wmn = 0. Meanwhile, the Bianchi 
identity V \lRpq)mn = leads to V m (Rmn — \9mnR + ^9mn) = 0. Thus, assuming that the 
parameter A and are constants, we have 

k a V m T 9 mn = \V M C MN . (21) 



1 If appropriate solutions exist in this system, we do not require the cosmological constant to be negative. 
2 It was also named as C-Tensor in [5]. 
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On the other hand from Eqs.(16) and Eqs.(19) one has 

V M T MN = v N (v 2 9 - ^) , V M C MN = V fRR, (22) 

from which together with Eqs.(21), we arrive at either — Vjv$ = or Wg = 0. Thus if we 
assume that the pseudo scalar field 9 is space-time coordinate dependent, i.e. vn ^ 0, the Bianchi 
identity together with Einstein equations will lead to the pseudo scalar equation W$ = 0. 



3.2 Non-relativistic Hydrodynamics Expansion 

From the above discussion we can see that the pseudo scalar equation is not independent. Thus, 
we will focus on the gravity equations henceforth. The Ricci scalar can be obtained from the trace 
of Eqs.(14) as 

R - 4A + XC% - K A [{d6) 2 + IV {9)} = 0, (23) 

which lead to the trace- reversed form of the gravity equations: Emn + ^cmn = 0, where 

E MN = Rmn — J^Qmn — K 4 t MN , t M N = d M 9d N 6 + gMNV(9). (24) 

The simplified Cotton tensor cmn = Cmn is due to its traceless property C% — 0. Then the 
equations of motion (EOMs) to be solved are 3 

Wmn = 0, Wmn = E M n + ACmtv, (25) 

We = 0, ^^fi-^-jM (26) 

do 4 

We assume that the following general (3+l)-dimensional black brane metric 

ds 2 = -g tt {r)dt 2 + g rr (r)dr 2 + g xx (r)5 ij dx l dx 3 , z,j = l,2, (27) 

as well as a radial coordinate dependent pseudo scalar field 9(r) solves the EOMs. With Dirichlet 
boundary condition at the finite cutoff surface r = r c , keeping the induced metric flat, and employing 
the non-relativistic hydrodynamics expansion method (see Appendix A for details), one can perturb 
the background metric up to the second order of expansion parameter e<l, 

d s !h\ = gniNdx M dx N = 2g r Jr)drdf — ^ / dr 2 + ^—j \ 5ijdx % dx 3 
(> gtt{r c ) gxx(r c ) 

- 2g r ,(r)(3 l dx*dr + 2 ( - p. d£ i df 

+ g r ,(r)/3 2 drdf - (^f\ - (^df 2 + ft/^xW) 

\gtt\ r c) gxx\ r c) J 

9rr(r),^ , , x ~ fg'tt(r) 9u(r c )\ ( 9tt{r) „ 2 



(Sf) g r r(r) rr drdf + (Sr) — — — # rf (r)drdf r~\ d ^ 

grr{r) \9tt{r) gtt{r c )J \ 9tt\r t 

9L(r) g'xx(r c )\ gxx{r) s JAijjM - , , n ^ 3 , 



(5r) _ vjcxkcj ^j_ 6lj dx*dx> + 0(e 3 ), (21 

\gxx(r) gxx{r c )J gxx(r c ) 



5 In what follows, we will use the normalization 2k 4 = 1, which means 16^6"^ = 1. 
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where we have defined g r f(r) = \J g rr (r ) -y/ gu (r) / gu (r c ) ■ The bulk coordinates are x M = (r, x M ), 
the intrinsic coordinates x M = (f , x l ) = (f , x, y) and we denote <9j = d x i A . The pseudo scalar field 
should also be expanded up to the second order 



0(r) = 9(f) = 9(r) + 9'(r)Sf. 



(29) 



We have assumed that the perturbation parameters are intrinsic coordinates x^ dependent, where 
Sf = f — r = 5f(x^) is a scalar perturbation at order e 2 and /3* = is the Lorentz boost 

parameter at order e 1 . Together with non-relativistic low frequency, long wavelength limit, we use 
the following scaling 



Or 



(30) 



Then we can solve the EOMs order by order with the non-relativistic hydrodynamics expansion. 



The Cotton tensor firstly presents at order e 2 , namely C 
background equations 



(o) 

MN 



and C 



(o) 



r(0) 



(1) 
MN 



0. Thus we obtain 



(31) 



where three of them are independent, and 
EB = 



E 



(o) 



2 \g rr {r 
9tt(r 



9tt{r) 



9:< 



.9., 



(0) 



gxx 



A 



f Iwfl(r)) + 1 ^ 

2 2g„(r) g tt (r) [g' tt (r) 

1 1 ^(r) 

2 2g rr (r) g xx (r) 



gxx 

9' xx ( 



[r) 
r) 



g* 



9'rr( r ) , 



9'L(r) 



9xx( r ) 



1 

" 2 

9' rr ( r ) _ 9u( r ) 

g rr (r) gtt(r) 



w, 



(0) 



0"( 



^LO) , g'tti r ) 



9rA r ) 



9* 



2gtt(r) 2g r 



9'(r) 

g rr (r) 



dV(0(r)) 
d9{r) 



(32) 
(33) 
(34) 
(35) 



We assume that background equations are always satisfied and they can be used to simplify the 
perturbation equations at higher order. It turns out that the intrinsic coordinate dependent metric 
and pseudo scalar can only solve the EOMs (25) and (26) up to order e 1 . Employing the spatial 
rotational 50(2) symmetry between x l directions, the nonzero traceless tensor sector and scalar 
sector appear in Wmn firstly at order e 2 . Correction terms to the metric and pseudo scalar are 
required to cancel the corresponding structure. For the traceless shear tensor sector, the following 
tensor correction terms are needed in the metric at order e 2 



9xx{r)y/9u(r c ) 



(2F(r) a tJ + 2F A (r) a*) dx l dx\ 



a 



'.i 



(36) 
(37) 



where F(r) and Fa{t) are chosen to cancel the ordinary shear tensor and Hall tensor af- of source 
terms at order e 2 , respectively. For the scalar sector, we will choose the gauge that no corrections 



'Here and henceforth, subscript i denote x or y. 
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are needed in the pseudo scalar, and we introduce the following terms in the metric to cancel the 
residual curl scalar ft appearing in the Cotton tensor at order e 2 



2g r f(r)H r f(r)drdf 



gtt[r) -H^{r)df 2 + ^P-H ££ (r)dx t dx l ) ft, 



9tt{r c 



Qxx (j"c) 



(38) 
(39) 



In summary, to solve the trace- reversed form of the gravity equations Wm n = and the pseudo 
scalar equation W$ = 0, we will work with the perturbed metric and pseudo scalar field with 
corrections up to e 2 in the non-relativistic hydrodynamic limit, 

ds 2 = ds\ b) +ds 2 {tye2) +dsf Sj£ 2 ) , (40) 

0(r) = 6(r) + 9'(r)Sr, (41) 

where the subscripts (b), (t), (s) in the metric represent "background", "tensor", and "scalar", re- 
spectively. 6{r) is substituted for 9 in EOMs. Correction terms to the vector sector appear at order 
e 3 , thus will not be considered in this paper. 



3.2.1 Traceless Tensor Sector of the Perturbation Equations 

By substituting the metric (40) and pseudo scalar field (41) into the EOMs Wmn 
tensor sector gives the following second order ordinary differential equation 



0, the traceless 



= — 



d_ 

dr 



+ — 



d_ 

dr 



9, 



9xx{r) 



V9tt(r) 
\j9rr(r) 



F'(r) + 1 



a 



'j 



F'Jr) + 



A 



9'u(r) 
2grr(r) \9tt(r) 




(42) 



As <7ij and have different tensor structure, we can solve the two second order differential equa- 
tions, respectively. Here we are interested in the physics between the black brane horizon and 
the cutoff surface S c . At the cutoff surface E c , we impose the Dirichlet boundary condition to keep 
the induced metric ^^ v {r c ) flat, i.e., F{r c ) = and Fa{t c ) = 0. At the horizon r = rn, we demand 
Fhh) and F^Orh) to be finite. Then the function F(r) and Fa(t) can be determined as 



F(r) 



F A (r) 



9tt(y) 

'9rr(y) 



9xx(y)\ ' 
A fg'uiv) 



9tt(y) [2g r r(y) \gu(y) 

where the two integration constants are 

Cf = 9xx{rh), 

X gxx(r h ) ( g' tt {r h ) g' xx (r h ) \ 



9xx(y) 
gxx(y) 



c Fa 



9xx(y) 



^gxx{r h )g'tt{ r h\ e ,^^ 



(43) 
(44) 

(45) 
(46) 



2 g rr {rh) \gtt(r h ) gxx(r h )J 2 g rr (r h ) g tt (r h ) 

Here we have used the assumption that at the horizon of the black brane solution (27), gu{f) has 
the first order zero gui^h) = 0, and g rT {r) has the first order pole g^i. r h) = 0, while their product 
gtt(r h )g rr (r h ) is finite. 
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3.2.2 Scalar Sector of the Perturbation Equations 

After solving the traceless tensor perturbation equations, only scalar sectors are left at order e 2 . 
At first, it is worthy to emphasize that the constraint equations on the scalar sector of gravity 
equations at order e 2 give 

'9'ttij) 



N M W Mf 



0. 



w (r) 



\9tt(r 



9* 



(47) 



where N is the unite normal vector of the constant r surface. As Wq(t) is not identically equal to 
zero, the constraint equation (47) leads to the incompressible condition G = dif3 l = 0, which shows 
the velocity field is divergence-free while the vorticity Q = e^d^j is still allowed. 
Then we need to solve scalar sector of the perturbation equations in (25) 



= 0, m 



(2) 



0. 



+ w 

xx 1 yy 



(2) 



0, w, 



(2) 



0. 



(48) 



Among them, only three are independent. Using the solutions (40) and (41), we can obtain the 
scalar sector of gravity equations defined in (24) at order e 2 , 



rr 



■HUr) + 



+ 



g'u(r) H' xx (r) 
g tt (r c ) 2g rr (r) 
9tt(r) g'^jr) 1 g' xx (r) 
9tt(r c ) g xx {r) 2g 



9rr( r ) , 9tt( r ) 

2g rr (r) 
HUr) 



2g tt (r) 
9tt(r) 
9tt(r c ) 
9u( r ) 



9a 
9xx 



HUr) - Hi 



xA r ) 



n. 



lg' rr (r)\ , 3g' tt (r)\Hi 



9xx{r) 2g„{r)) 2 g tt (r c ) J 2g rr (r) 
9u(r) g' rr (r)\ \ HMr)-2HMr) 



9tt{r) H'; T {r) 
9tt(r c ) 2g rr {r) 



E 



(2) 



E 



(2) 



yy 
9xx(r) 



J XX 



XX yr) \gtt{r c ) gtt(r c ) 9rr y 

U r ) H 'rA r ) ~ H 'rr( r ) , ( 1 9xx(r) ( g' rr ( r ) 



<xx(r c ) 9rr(r) 

H'L(r) , ( sf„{r) fg'M 



4g X x(r c 
_ 9'tt( r ) 

9xx(r c ) 9rr(r) ' \2g xx {r c ) \g rr {r) g u {r) 



2g rr (r) 
9'tt{ r )\ 9'xx( r ) 



n. 



9rr(r) g tt (r) 



9., 



9l(r) 



H, 



2H r 



9rr{r) 



9rr{r) 

n. 



(49) 



and we have imposed the incompressible condition here and henceforth. The corresponding scalar 
sector of Cotton tensor defined in Eq.(17), appears firstly at order e 2 , can be expressed as 



C (2) 



1 



9tt(r) g X x(r c ) g rr (r) 



C 



(2) 



'g rr [r) V 9tt(r c ) 9xx 
-1 9tt{r) g xx (r c ) d_ 
2^v 7 )gT{r c )gl / x\r) dr 



d_ 
dr 



gU r ) 

g rr (r) \g a (r) 



6'(r) 



n. 



^gtt(r)^g xx (r) f g' tt {r 



r (2) r (2) 
xx ' yy 



grr 
d 



2^g rr {r) y/gu(r c ) y/gu(r) dr 



T) \9tt(r) 
gtt{r) ( g' tt {r) 



9'(r) 



n. 



.9rr(r) \gtt(r 



gL( r ) 



e'(r) 



n. 



In the end, we have the pseudo scalar equation defined in Eq.(26) at order e 2 , 



W, 



(2) 



fl'(r) 

_g rr (r) 

A 2 



(HL(r) - H' r , 



H' xx (r)) + (HMr) ~ 2H rf (r)) 



d0{r) 



gxx(r c ) g' xx (r) g tt (r) ( g' tt (r) g' xx (r) 



4 gU 2 (r) gxxir) g xx {r) y g u (r c ) \g u (r 



n. 



(50) 



(51) 



With Dirichlet boundary condition at the cutoff surface, namely Hf?(r c ) = 0, H xx (r c ) = and 
the finite boundary condition at the horizon, the function H^(r), H r ^(r) and H^(r) are expected 
to be solved. Similarly procedure could be carried out by following the approach in Ref. [26], here 
we just assume these formal functions solving the scalar perturbation equations (48), similar to the 
assumption the black brane metric (27) solving the background equations. 



3.3 Incompressible Navier-Stokes at the Cutoff Surface 

We now study the dual (2+l)-dimensional hydrodynamics in non-relativistic limit to the background 
gravity solution. At the cutoff surface E c , the Brown- York tensor [27] 

% u = 2 (K% v (r c ) - + Cv(r c )) (52) 

is identified with the stress energy of the dual fluid in Ref. [32], where is the extrinsic curvature 
tensor, ^^ v {r c ) = r]^ is the induced metric on the cutoff surface, and C is a constant. In the 
non-relativistic limit, the non-zero components at order e° turn out to be 

7#> = e (r c ), e (r c ) = -^l= 9 -^\-2C, (53) 

\/9rr{r c ) 9xx{r c ) 

if = Po(r c )^, Po (r c ) = -J= + ^4) + 2C. (54) 

yjg„{r c ) \9tt{r c ) 9xx{r c ) J 

eo(r c ) and po(r c ) could be regarded as the energy density and pressure of the dual fluid at rest 
on the cutoff surface in the thermal equilibrium state. They are related with each other through 
the equation of state of the fluid, which comes from the Hamiltonian constraint (see Eq.(106) in 
Appendix A). At the zeroth order it is given by 

1 (e (r c ) + 2C) (e (r c ) + 4p (r c ) - 6 C) = A+ M + i Q - —Lj) 6'\r c ) (55) 

The entropy density and temperature of the dual fluid is defined as the local physical quantities at 
the cutoff surface S c of the black brane [32, 39] 

So(r c ) = ^^44 T(r c ) = T H = lim - /f } (56) 

4Giv gxx{r c ) V9tt{r c ) ^r h Ait^/ g tt {r)g rr {r) 

where a factor in the entropy density so(r c ) has been chosen to meet with the Bekenstein-Hawking 
entropy formula, while the local temperature at the cutoff surface T(r c ) relates with the Hawking 
temperature T# through the Tolman redshift relation. 

At order e 1 , only vector components appear in the dual stress energy tensor 

Tff = 7-« = -Mrjfk, Mr c ) = e (r c ) + Po (r c ) = -J== ~ ^) , (57) 

y/g r r(r c ) \9tt(r c ) g xx {r c ) J 

where Wo(r c ) is the heat function per unit volume. While at order e 2 , the scalar and tensor modes 
are 

if = w (r c ) (S + (3 2 ) + (' A (r c )Q, (58) 
%f = w (r c )(V + f3 l f3 J )-2r l (r c )a l] -2 VA (r c )o? 3 -CA(r c )tt5 lv (59) 
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The normalized energy density and pressure perturbations are denned as 



_ de {r c ) 5r{x^) _ dp (r c ) 5f(x^) 

dr c wo^c)' dr c w (r c ) ' 

where 5f(x fl ) could be regarded as the thermodynamic parameter due to the temperature pertur- 
bation. 

The transport coefficient appears at the first order gradient expansion. In the traceless tensor 
modes, the normal shear viscosity reads 




I 9tt{r e ) . g xx (r h ) 7](r c ) 1 

Wc) = l + \ — r ^F{r c ) = — , —-=4G N = — , (61) 

grr{r c ) 9xx\J c ) s {r c ) An 

which shows that tj/sq is cutoff scale independent. While the Hall viscosity turns out to be 



, ytt{r c ) , _ A g xx (r h ) 6'(r h ) g' a (r h ) A w (r c )9'(r c ) 
J 9rr(r c ) A c) 2 g xx (r c ) g rr (r h ) g tt (r h ) 2 rf^j > ^ 



Va(t c ) A 9'{r h ) g' tt (r h ) A g xx (r c ) 6'{r c ) f g' tt (r c ) g' xx (r c ) 



s (r c ) 8ir g rr (r h ) g tt (r h ) 8ir g xx (r h ) g rr (r c ) \g u (r c ) g xx (r, 



(63) 



Thus, we see that t/a/so is cutoff scale dependent due to the second term. When the curoff surface 
approaches to the horizon, namely, r c — > r h , one has f]A(r c ) / s (r c ) — > 0. Even without solving 
the scalar perturbation equations in (48), we can still give the formula of curl viscosity with the 
formal functions in Eq.(38). In the dual fluid stress energy tensor (58) and (59), the scalar sector 
coefficients at order e 2 are 

C A (r c ) = f^^firf (r c ) - HUrS) , (64) 

Ca(t c ) = -=L= ( 9 -^H rf (r c ) - HL(r c )) + (65) 
y/9rr{r c ) \9tt{r c ) J 2 

Here we have used the Dirichlet boundary condition H?f(r c ) = and H xx (r c ) = 0. If the Landau 
frame is chosen that C/i( r c) = 0, we can identify Ca(^c) as the curl viscosity of the dual fluid. 

According to the Gauss-Codazi equation in pure geometry sector (see Eq.(107) in Appendix. A) 
and the gravity equations in Eqs.(14), we have 

d»%u = -2N A R AB h B u = N A (2\C AB - T e AB ) h B u , (66) 

where h^s = g A B — N A N B is the induced metric at the cutoff surface. In our gauge choice of the 
scalar perturbations in Eq.(40)-(41), it could be checked that up to order e 3 of the source terms on 
the right-handed side, the Cotton tensor Cab contributes nothing, while the stress energy tensor 
T 6 AB contributes an external force density as 

V9rr{r c ) 
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which comes from the scaling transformation parameter 8f[x^), and the normalized force density 
is fi( r c) = Fi{ r c) /wq{t c ). Thus finally we get the incompressible Navier-Stokes equations with Hall 
and curl viscosities in Eq.(lO) from Einstein gravity with Chern-Simons term, 

dS % = 0, d f fr + fr&fr + d(P - v{r c ) d 2 fr - v A (r c yid 2 (3 3 - U^d^^ = f z (r c ). (68) 

where the kinematic viscosities are defined as 

K r c) = rj(r c )/wo{r c ), v A {r c ) = VA(r c ) /w {r c ), £a(»"c) = (A{r c )/w (r c ), (69) 

respectively. The Reynolds number of the dual fluid associated with the shear viscosity and Hall 
viscosity can also be defined as 

uL 1 uL 1 

K = OC —f-^r, K-A = OC — , (70) 

v u(r c ) v A v A [r c ) 

where u is the characteristic velocity, L is the characteristic scale. If we further set 5r(x fJ ') = 0, 
then the normalized pressure perturbation V and force density fi(r c ) vanishes, the holographic 
incompressible Navier-Stokes equations will reduce into the incompressible Burger's equations [28] 
with Hall viscosity and curl viscosity, 

d t (3 l = 0, dffii + fatf - v{r c ) d 2 (3 t - v A {r c )e ij ^^ - U{r c )e jk d t d 3 (3 k = 0. (71) 



4 Dual Fluid at the Infinite Boundary 

To compare with results in Ref. [5], we take the limit with the cutoff surface r c — > oo. The ansatz 
of the bulk metric is 

ds 2 = -r 2 f(r)dr 2 + 2H{r)drdr + r 2 dx i dx i , i,j = 1,2, (72) 

and 9 = 9{r) is a pseudo scalar coupled to the gravitational Chern-Simons term, as well as a 
negative cosmological constant A = —3/£ 2 . Under this ansatz, the Pontryagin density RR is 
identically zero [11]. This metric is related to our generic metric in Eq.(87) through 

H(r) 2 

gtt(r) = r 2 f(r), g xx (r) = r 2 , g rr (r) = (73) 

r 2 j{r) 

And the pseudo scalar equation Wq = 0, with Eq.(35) becomes 

K ' \r fir) H(r)J K ) r 2 f(r) d9(r) K ' 

The dynamics shear viscosity and Hall viscosity at the cutoff surface S c defined in Eq.(61) and (62) 
turn out to be 

, x 1 Th ( v A 1 f rjf(r h )e'(r h ) rjf{r c )e'{r c ) \ 



lfaG N r c ' " v c> lQnG N r 2 c \ 2H(r h ) 2 2H(r c ) 2 
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According to Eqs.(56), the entropy density and temperature of the dual fluid at cutoff surface E c 



are 



i \ 1 Th t/ ^ r hf'( r h) ,„„•, 

s (r c ) = — , T(r c ) = ==, T H = ■ 76 

AG N r c r c y/f(r c ) A-KH{r h ) 



Thus the shear and Hall viscosity over entropy density are given respectively by 
r){r c ) 1 VA(r c ) A 1 (rif{r h )6\r h ) rif'(r c )e'(r c ) 



s (r c ) 4vr ' s {r c ) 4n rf t \ 2H{r h y 2#(r c ) 2 



(77) 



where we have used normalization 167^^ = 1- The shear viscosity over entropy density of dual 
fluid does not run with the cutoff surface, while it is not the case for the Hall viscosity. In the 
infinity boundary limit r c — > oo, if we make the following assumptions that 



f{r c ) -+ 1 - 0{r 6 h /r% H(r c ) -+ 1 - 0(r£/r c d ), 6(r c ) ^ 0(r»/r c «), (78) 

when a > 0, we can drop the second term in tja/s, thus recover the result in Ref.[5] 5 . In that 
case, the Hall viscosity over entropy density is entirely determined by the near horizon region of 
the black-brane geometry. 



5 Discussions 

In this paper, we have investigated the non-relativistic hydrodynamics with Hall viscosity and curl 
viscosity living on a finite cutoff surface, dual to (3 + l)-dimensional Einstein gravity with a pseudo 
scalar coupled to a topological gravitational Chern-Simmons term. The topological Pontryagin 
density RR is totally depicted by the Weyl tensor, which describers the traceless part of the Riemann 
tensor. While the Ricci tensor is Riemann tensor's trace and together with Weyl tensor provide a 
complete description of the curvature of the space-time. In (3+l)-dimensional space-time, the Weyl 
tensor and its dual can be written as 

Cabcd = Rabcd — {9a[cRd\b — 9a[dRc]b) + -^9a[c9d]bR, (79) 

riABCD 1 CDEF n AB fin — fiABCDri /orA 

G = 2 EF ' ~ ^BACD, {OVj 

and it could be shown that the Pontryagin density RR = CC (for example see [13]). One can 
further define the gravito-electric and gravito-magnetic field as 

E AC + iB AC = {Cabcd + iC A bcd)u b u d = {Cabcd + iC A bcd)P BD , (81) 

where u B is a normalized time-like 4-velocity, and its projection P BD = g BD + u B u D . The last 
identity is true since the traceless properties of the Wely tensor. Thus, the Lagrangian density in 
Eq.(12) in the bulk can be expressed as 

C cs = j ORR = ^6CC = -AX 6 E AC B AC , (82) 

J Thcrc is a notation difference from Ref. [5]: Ai lcrc = — Athoro- 
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which is analogous to that in U(l) gauge theory, Cq S = kOF^F^ 



- -4k6(E ■ B). With suitable 
coefficient and a constant ^-vacuum, Cq S can be used to describe the electromagnetic response of 
(3+l)-dimensional topological insulators (TI), while Ccs can be used to describe the gravitational 
response of (3+l)-dimensional topological superconductors (TSC), a half quantized thermal Hall 
effect may appear on the surface [23, 24]. Alternatively, if we consider the stress response of 
(3+l)-dimensional TI to an external torsion field, it was studied in Ref.[25] that the time-reversal 
invariant TI will exhibit a quantum Hall viscosity on their surfaces. Thus it would be interesting to 
see whether the bulk gravity in this paper could describe a deformed TI or TSC, with the surface 
quantum Hall effect as the dual description. In our general perturbative metric (28), the Lagrangian 
of the anomaly term firstly appears at order e 2 as 



r X ar„\ 9xx(r c ) g' xx (r) / g u {r\ 




S3) 



which has been simplified through the background equations in (31). To get a non- vanishing 
Hall viscosity t)a in Eq.(62), the pseudo scalar field 9{r) is required to be coordinate r-dependent. 
Otherwise, when 9{r) vanishes, parity is not broken in the bulk and in the dual boundary theory, 
and when 9(r) is a non-zero constant, the gravitational Chern-Simons term is just a surface term in 
the action which violates the parity in the bulk but contributes nothing to the equations of motion 
in the bulk, hence leads to a vanishing Hall viscosity. In Ref.[6], by including a gauge field in the 
bulk gravity such that the solution is free from the violation of the positive energy theorem [29], 
the pseudo scalar hair would break parity spontaneously. 

Inheriting the spirit of holographic Wilson's renormalization group (RG) approach in Ref. [32], 
the re-scaled cutoff size could be regarded as the RG running scale. The non-vanishing trace of 
the symmetric stress energy tensor leads to the trace anomaly in quantum field theory, where the 
corresponding /3-function can be obtained from the coefficients in front of the trace. In our dual 
gravity model, Setting A = to turn off the gravitational Chern-Simons term in the bulk, we can 
go back to the d-dimensional fluid corresponding to (d+l)-dimensional bulk gravity with a scalar 
field 9(r), the relation between the trace in Eq.(105) of dual Brown- York stress energy tensor 
at zeroth order and /3( e ) is, 



"/3(e)- 



d 







d/2, 

9xx [r c 



9xx ( r c 



y/9rr{rc) 



/3(e) = r c — [e(r c )g d J x 2 (r c )] 



M) 



where e(r c ) is given in Eq.(103), r c is analogous to the cutoff energy scale. Considering there exists 
a conformal factor difference in the definition of the dual stress tensor in Ref. [37] and transforming 
to that frame, we can conclude that the trace of the fluid stress tensor is generated by the energy 
density as well as an extra term contributed by the the scalar field. This relation could also be 
generalized at higher order of hydrodynamic expansion, and in order to consider the effects of parity 
violating interactions, one can turning on the gravitational Chern-Simons term in the bulk. As the 
ratio of the Hall viscosity over entropy density in Eq.(63) is found to be cutoff scale dependent, we 
can also define the following /3-function 



R 9 



VA{r c 



Xr c d 
8% dr r 



9xx(r c ) w (r c )9'(r c ) 
9xx (r h ) y/grr(r c ) 



55) 
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to represent its non-trivial evolution, where wo(r c ) is the heat function per unit volume in Eq.(57). 
When the cutoff surface approaches to the horizon of the background spacetime, rj{r h ) / SQ{r h ) = 0, 
but 0( VA / S ) could be non-zero. 
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A Non-relativistic Hydrodynamic from Gravity 

This appendix briefly introduces the non-relativistic expansion procedure associated with a cutoff 
surface [39]. For further details, see Refs. [30, 31, 32, 33, 34, 35, 36, 37, 38]. 

A.l Bulk Geometry and a Finite Cutoff Surface 

In order to study the fluid in d- dimensional flat space-time, let us consider a (d + l)-dimensional 
background geometry with a generic metric 

ds 2 d+1 = -g tt {r)dt 2 + g^ir^^dx'dx 3 + g rr (r)dr 2 , i,j = 1, . . . d - 1, (86) 

where the metric components depend only on radial coordinate r. In addition, we assume the 
fully spacial rotational SO(d — 1) symmetry in x l directions, and require that the geometry of the 
space-time manifold M. has a well-defined future horizon located at r = rh, where gtt{f) has the 
first-order zero gttijh) = 0, and g Tr {r) has the first-order pole grr( r h) = [30]. By using the 
Eddington-Finkelstein coordinate r defined through dr = dt+ \/ g rr ( r ) / 9tt( r ) dr, we can rewrite the 
bulk metric as 

ds 2 d+l = -g tt {r)dT 2 + g xx (r)dxidx l + 2^g tt (r)g rr (r)dTdr. (87) 

Let us define an arbitrary hyper-surface at a constant radial coordinate r as S r , and introduce 
a specific finite cutoff hyper-surface S c at r = r c outside the horizon (r c > r^), where the associated 
intrinsic coordinates x M ~ (f, x % ) on S c are 

x° = t = ry/g u (r c ), x l = iVfeW, i = l,2,...,d-l. (88) 
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The bulk metric in Eq.(87) in the intrinsic coordinates becomes 



ds 2 d+l = -^f\df 2 + ^P-dxdx 1 + 2 J ^^Jg~Jr)drdr. (89) 
9tt{r c ) 9xx(r c ) ]] 9tt{r c ) 

We adopt the Arnowitt-Deser-Misner(ADM) decomposition along r direction [40] 

ds 2 d+1 = a 2 dr 2 + ^(dx" - ^dr)(dx u - fi u dr), (90) 

with the "lapse function" a = g rr ( r ) and the "shift vector" /3 M = gtt( r c) / 9tt( r ) \J 9rr( r ) 
respectively. The generic induced metric 7 M i,(r) is an analytically extension from S c to £ r , where 

1 „ v (r)dx»dx v = -^fXdf 2 + ^Idx.dx 1 . (91) 
' 9u(r c ) 9xx{r c ) V ' 

It is worthy to notice that 7 At „(r c ) = r]^ = (—1,5^), which implies that the induced metric 7^(7") 
on E r reduces to a flat metric on £ c . 

A. 2 Diffeomorphisms Associated with the Cutoff Surface 

With a flat induced metric given at the cutoff surface £ c , we can take two finite diffeomorphisms. 
The first one is a linear scale transformation along the radial coordinate r and the rescaling of the 
intrinsic coordinates = (f,x l ), 



r^r = k(r), r^rJ 9 -^ x^ x\ 9 -^\, (92) 

V 9tt{r c ) V 9xx{r c ) 

where k(r) is a linear function of r, and a concrete form of k(r) will be chosen according to the 
specific global geometry of the bulk metric in Eq.(89). In addition, the notation r c = k(r c ) is 
introduced. Under this diffeomorphism, the bulk metric in Eq.(89) is rescaled as 



9u^) d ^ + 9^(r) d£id ^ + 2 9uif) 
9tt{r c ) 9xx{r c ) y 

The second diffeomorphism is a generic Lorentz transformation 



ds 2 ds 2 = -^-df 2 + ^J-dxM 1 + 2J^{y/^(fjdfdf. (93) 
9tt{r c ) 9xx{r c ) V 9tt{r c ) 



r 7 (f - # # - 10 + (7 - i)^r^', (94) 



with the boost parameter (3j = 5ij(3 l and Lorentz factor 7 = — /3 2 . We can define the 

c?-dimensional velocity, and <i- dimensional polarization vector as 

% = 7(-l, A), n> = (-7^, ^ + (7 - 1)^), (95) 

with u^Uu = —1, n^nf, = d — 1 and u M n* = 0. Then the Lorentz transformation on the cutoff 
hypersurface becomes f — > — and — >■ n^a; M . The re-scaled bulk metric in Eq.(93) becomes 



ds 2 d+l ->• ds 2 = -^^-u.u.dx^dr + ^P-P, u dx^dx u - 2 J ^^^Mu.dx^dr, (96) 
9tt{r c ) 9xx{r c ) V #4 r c) 
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where the projection tensor 

Pfiu = f)iiv + U^Uv = SijU^n^ = ^ _^,2^i fii _|_ jZflip . ^ ' (^7) 

satisfying — d — 1. Under the two diffeomorphism transformations, the metric in the ADM-like 
decomposition with the "lapse" function a(r) and the "shift" vector (3(r) becomes 

ds 2 d+l = a 2 {r)df 2 + % v (r) (dx» - 0»(r)dr) [dx v - $"(r)drj , (98) 



V 9tt{r) 

and the induced metric in Eq.(91) is transformed into j^r) as 



(99) 



%u(r) = -^u,u v + ^\P, U . (100) 
9u{r c ) 9xx\r c ) 

At the cutoff r = r c surface, which is equivalent to f = f c , the induced metric 7 M ^(r c ) keeps flat. 



A. 3 Dual Hydrodynamic in Non-Relativist ic Limit 

We perturb the diffeomorphic metric (96) in the non-relativistic limit. As r is a radial coordinate 
along the extra dimension orthogonal to the cutoff hyper- surface, we can define a scaling transfor- 
mation parameter 5r = r — r. Assuming the scalar parameter Sf and the velocity parameters /?* in 
the Lorentz boost are all functions of the intrinsic coordinates x^ = (f,x l ), in the non-relativistic 
hydrodynamic limit we have the following scaling ( here the subscript i present x 1 ) 

d r ~ e°, di ~ ~ e 1 , d f ~ 5r{x>") ~ e 2 , (101) 

where e C 1 is a small parameter. Therefore, the metric function can be expanded, i.e., gttif) = 
9tt{f) + 9tt{ r )^^{ r ) + • • •> then the diffeomorphic metric could be expanded order by order in the 
expansion parameter, up to order e 2 , just gives the result in Eq.(28). In Einstein gravity, we assume 
that the generic metric in (86) solves Einstein's equations. Then the diffeomorphic metric is also the 
solution. After demanding the transformation parameters and 8f to be coordinates x^ dependent, 
the metric in (96) is no longer the solution of Einstein's equations. We can solve the Einstein's 
equations order by order in the non-relativistic hydrodynamic expansion, via adding new correction 
terms to the metric. Similar procedure could be carried out when the matter field are present. 

The intrinsic curvature of the bulk geometry is measured by Riemann tensor, while the extrinsic 
curvature depends on how the hyper-surface is embedded into the bulk geometry. Thus, when we 
perturb the background metric, we can choose the gauge to keep the induce metric ^^{r^} flat. The 
extrinsic curvature at the cutoff surface is K^ v = \Cf]^$(r)\ r=rc , where Cjy is the Lie derivative 

along the outpointing normal vector N of the cutoff surface, 7^(r) is an analytic extension of the 
induced metric 7 M! y( r c) on E c . And they are all associated with the perturbed metric with correction 
terms. In (d+ l)-dimensional Einstein gravity, if we use the normalization condition IQttG^ 1 ^ = 1, 
the Brown- York tensor on the cutoff surface S c is given by [27] 

Tf = 2(K% u (r c ) - K„„ + CvW). (102) 
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It is identified with the stress tensor of the dual fluid in Ref. [32]. In the non-relativistic hydrody- 
namic expansion, we can also fix the physical stress energy tensor order by order. The zeroth order 
stress tensor could be written as 7^ y = eir^u^u 1 " + p(r c )P^ u , where 

e(r c ) = -J^ 9 -^ 2C , (103) 

pW = 7m(^ +id - 2) iM) +2C < (104) 

are the energy density and pressure of the dual fluid, respectively. And the trace of the dual stress 
energy tensor at the zeroth order is 

T£ = -e(r c ) + (d-l)p(r c ). (105) 

In addition, the dual fluid is described by the constraint equations [41, 32] 

(k 2 -K AB K AB ) = 2G MN N M N N \^ (106) 
2D A (kh AB -K AB ) = -2G MN N M h N B \ Xc . (107) 

where Gmn is the Einstein tensor, h AB = g AB — N A N B and D denotes the derivative operator 
associated with the induced metric h AB . The first one is the Hamiltonian constraint which gives the 
equation of state of fluid relating the pressure and energy density. The second one is the momentum 
constraint which gives the evolution equations of the fluid, and it reduces to the incompressible 
Navier-Stokes equations in the non-relativistic limit. Correction terms will appear if matters or 
higher order derivative terms are added to the Einstein gravity. 
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